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SUMMARY 


The electronic polarizability « has been determined by perturbation methods for Li, ret a 
Ca and Se. The numerical calculations have been performed on the Swedish electronic digital 
computor BESK. The calculated values of « (« in units 10-"* em®) are: 20 for Li, 0.83 for Cat 
52-153 for Ca and 154 for Sc. Although these values of « probably are not very exact, they 
are however certainly of the correct order of magnitude. 


Introduction 


A great many experimental investigations have been performed to determine 
the electronic polarizability « of ions [1]. In most cases, « has been determined 
for ions with rare-gas-like configurations. However, Drechsler, Miiller and Henkel 
[2], [3] have worked out a method for measuring « of other species with the 
field emission microscope. They have published experimental values of « for Li, 
Ba, [2] and some metallic atoms [4]. 

Kirkwood’s variation method has been used several times for theoretical de- 
terminations of «. This is the case with works by e.g. Buckingham [5] and Guy 
and Harrand [6]. Pople and Schofield [7] have recently improved the method 
and applied it to argon. Moreover perturbation methods have been used by Stern- 
heimer [1], [8], who has calculated « for atoms with rare-gas-like configurations, 
i.e. s2p% and s? p$d1°, 

In a preliminary report [9] we have pointed out that it would be of great 
interest to find even rather approximate values for the electronic polarizability 
a of atoms and ions with partly filled d-shells. Therefore an attempt was made 
to find the order of magnitude of this polarizability from a simple perturbation 
calculation. Atoms with one, two or three electrons outside the rare-gas-like con- 
figurations were considered and the result was applied to lithium 1 s?2s, calcium 
4s? and scandium 3d4s8°. 


Theory 
The Hamiltonian of an atom in a uniform field F ist 
H=H,+ En cos 3; (1) 
1 Atomic units are used in the formulae throughout the paper. 
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where 7; and ®% are spherical polar coordinates of the i™ electron. The polar 
axis has the direction of the electric field. H, is the Hamiltonian of the free 
atom. 

When treating the effect of the electric field, it is first necessary to decide 


how the order of magnitude for the perturbing potential, # > 7; cos 0;, compares 
v 


with the other terms in the Hamiltonian [10]. Here the case is considered where 
the effect of the electric field is large compared with the spin-orbit interactions 
but small compared with the Coulomb interactions. As unperturbed Hamiltonian 
we then have: 


He=3(-4A.-Z) 452 (2) 


G<i Vij 


With this Hamiltonian the energy levels are classified by terms represented 
by the symbol ?°*Z. In the presence of the uniform field H we have axial 
symmetry and each term ?°*1Z splits up, so that states with different | MZ ,| have 
different energy values. In the case of scandium, the ground state is 3d4s?2D, 
The M, values are directly given by the m; values of the 3d-electron: (m= 
=0,+1,+2). The wave function for a special value of M, can thus be written 
as a single Slater determinant of unperturbed one-electron wave functions, 9; (Tr: €¢), 
where r; is the space coordinate and ¢; is the spin coordinate of the t™ elec- 
tron. We can thus consider our problem as a perturbed Hartree-Fock problem. 

The problem is now to determine the first-order perturbation gj of the func- 
tions y;. The energy can be determined up to the second order from q; and qj. 
It can be shown [11] that the first and second order perturbation energies take 
the form: 


E,=> [pf Er cos 8pidr=0 (3) 
E,- XJ gt Er cos Ig; dt (4) 
where g; is given by the equation: 
(—14.-2- a) + S0Glalie)—Glal “1 - 
= — Bry cos Fags (2)— 2G" lal) + Glal7’a)—-G'laléf)— Gl alis +e gi(2) 6) 
In 6) alii) = foi") gs (hy eu(2)dm ete 


&i is the orbital energy. By multiplying equation (5) by gj and integrating, 
€ is given by: 


= DUG tl gli)+ Gilal7i)—G' él glig)— Gélali7’y] (6) 
where (i| qj) = J 95" () gf (2) — 9, (I) qu(2)dr, dr, ete. 
12 
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The fact that €,=0 is equivalent to the fact that the atom has no permanent 
electric dipole moment [12]. As is well known, the electronic polarizability « is 
defined by the equation 


€,= —ta (7) 


The induced electric dipole moment is « Z. 

The purpose of this investigation was to find the order of magnitude of «. 
As it would be very laborious to work with the equations (5), the author has 
used an approximation of (5), namely: 


(-—4A+V,.—«)9i= —Er cos 39; (8) 
where 
(—4A+V,—«)gi:=0 (9) 


For given functions y;, V, is given by (9). In (8), both the sum and the last 
term on the right hand side of (5) are neglected, i.e. the influence on the per- 
turbation of gy; from the perturbations of the other orbitals is neglected. This 
approximation will be discussed further in connection with the discussion of the 
results. 

In the Hartree-Fock approximation the functions g; are of the form: 


a . 
P= Prin = YP (8, 9) (0). (10) 


From equation (8) and from the properties of Y7 (#,) [13] it is seen that g; will 
be of the form 


gi = (Risa (1) Via (8, p) + Bra (r) Vi" (8, p)) f (6). (11) 
Ri.1(r) and R,_,(r) can be written: 
is. (d+1—m)(l+1+m). 
Rizr) =— Cy] (ee co (12) 
: Rae (L—m) (I+ m) 
Ry-1(r) = 1 Pon) me (13) 


/ . . 
where P;,;,;, and P;_,;-; are given by the equations 


d? 1+1)(1+2 ' 
| fat 27 Bai] Pinta Pas (14) 
d2 l—1)l , 
|-aa ~ +2V)— 26] Pa=rPan (15) 


as has been shown by Sternheimer [1]. 
In the case of the scandium atom one must determine which of the | m|-values 
assigned to the 3d-orbital gives the lowest energy level. The contribution to é; 


from 3d for different |m| is: 
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|m | | E, (3d) 


where B= | Psa+s7Psadr (16) 
0 


and y= PsasptPsadr (17) 
0 


Now the calculation shows that # and y are both positive. This means that 
|m|=0 gives the lowest energy level. 


Numerical calculation 


The equations (14) and (15) have been solved by numerical methods using 
the Swedish electronic digital computer BESK. (For further details, see appendix.) 
As far as the author knows, there are no exact determinations of the ground 
state orbitals for Sc. As unperturbed one-electron orbitals, the author therefore 
has simply used orbitals of the analytic form r”e-“’, obtained from Slater’s rules. 
Thus the BESK-programme has been constructed for such orbitals. It is, how- 
ever, easy to modify the programme so that orbitals which are numerically given 
can be used. Although numerical SCF-orbitals are available in the case of Li [14], 
Ca** [15] and Ca [15], orbitals of the above-mentioned form have been used here 
also in these cases. The value of ~ was, however, obtained by fitting the ana- 
lytical orbital to the numerical one. In most cases, the exponent n was chosen 
to be an integer for computational reasons. In the case of the 4s-orbital of Ca, 
nm was also chosen equal to the Slater-rule n. 

As a check on the programme the polarizability of Ca**, previously calculated 
by Sternheimer [1], was computed. As is seen from Table 2, the agreement is 
satisfactory. 

The computed value of « can be expected to depend strongly on the choice of 
the unperturbed functions q; [16] [17]. For an estimation of this dependence the 
polarizability « of Ca 4s? was calculated for orbitals from Slater’s rules and for 
analytical orbitals chosen to fit the best SCF-orbitals. In Table 1, those values 


of « and n for Li, Ca** and Ca, which are chosen to fit the best SCF-orbitals, 
are given. 


Results 


Results for the polarizabilities of Li, Ca**, Ca and Se are given in Table 2. 
It is seen that « of Ca and « of Sc are both much larger than « of Ca**. 
In Table 3, the contributions to « from the various shells are given for Li 
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Table 1. Values of w and n for Li, Cat* and Ca. « chosen to fit 
the best SCF-orbitals. 


Atom Li Cat+ Ca 
Orbital a n yu n Le n 
ls 2.478 1 
28 0.667 2 
2p - — 
38 = ead .00 3 3.00 3 
3 p =e SS 2.83 3 2.83 3 
48 a — — 1.08 4 
48 2 a a cam 0.987 3.7 


Table 2. Calculated and experimental values of the electronic polarizability «. 


Atom | (10-2 em?) Method Author 
4 16+3 experimental Drechsler and Miiller [2] 
Li 20 theor., 4 chosen to fit the best Present investigation 
SCF orb. 
0.47-1.1 experimental Tessmann, see [1] 
Catt 0.73 theor., SCF orb. Sternheimer [1] 
0.83 theor., 4 chosen to fit the best Present investigation 
SCF orb. 
153 theor., 4 from Slater’s rules Present investigation 
Ca 52 theor., “ chosen to fit SCF orb. af 
57 ” ” 
Se 154* theor., 4 from Slater’s rules Present investigation 


* By mistake, incorrect figures were given for the values of « for Ca (Slater) and Sc in the 
preliminary report [9]. 


Table 3. The contribution to « from the various shells. « in 10~** cm3. 


\, Atom ’ 

(nl) Sc Li Se 
(1s) 0.03 10-5 
(28) 20 0.002 
(2p) — 0.003 
(33) = 0.24 
(3p) = 0.28 
(3d) = 26 
(48) — 127 
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Table 4. The contribution to « of Ca from the 4s shell for three different choices 
of the unperturbed 4s-orbital. « in 10% cm’. 


Ca 
Orbital n (48) fb (48) a (48) 
Slater’srule.... . 3.7 0.77 152 
Approx.of SCF ... 4 1.08 51 
Approx. of SCF 3.7 0.987 56 


and Sc. In the case of Li the contribution to « from the inner shell is negli- 
gible. In the case of Sc the same is true, ie. the main contribution originates 
from 3d and 438?. 

Table 4 shows the contribution to « of Ca from the 4s shell for different 
choices of the unperturbed 4s-orbital. 


Discussion 


From Table 2 it is seen that the author’s calculated value for « of Ca** lies 
within the range of the observed values. The value is slightly higher than Stern- 
heimer’s. The difference may be due to the fact that he used tabulated Hartree- 
Fock functions instead of which the present author used an analytical approx- 
imation. On the other hand, Sternheimer estimates the error in the values of « 
due to the interpolation procedure used to obtain gy; to be less than + 20%, 
whereas the present author estimates this error to be less than + 1% (see ap- 
pendix). The result for Ca** shows that the BESK-programme functions satis- 
factorily. 

From the results for Ca (Table 4) it is seen that the computed value of « is 
very sensitive to the choice of the unperturbed atomic orbitals. The value of « 
for the Slater’s rule orbital differs by a factor of 3 from the values for the two 
orbitals chosen to fit SCF-orbitals. 

It is difficult to estimate a numerical value for the error due to the fact that 
the approximate equation (8) has been used instead of the equation (5), but a 
rough estimate shows that this approximation would cause about the same per- 
centage error in « for all calculated cases. Now the calculated values of « for 
Ca** and Li are in good agreement with the experimental values, which could 
justify the approximation made. 

The result of this investigation supports the assumption that the polarizabilities 
of calcium and scandium are considerably larger than the polarizabilities of closed 
shell atoms (for which detailed theoretical investigations already exist). 

It is interesting to compare the polarizabilities of the s-electrons of different 
atoms. Some values are given in Table 5. It is seen that « increases with n and 
decreases with increasing charge within a certain period. It seems, however, as 


if the polarizability of H is unique, being much smaller than that for other atoms 
with valence electrons. 
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Table 5. Polarizabilities of one ns-electron. 


a 
Electron Atom 
Slater orb. SCF orb. 
ls H — 0.667 
28 Li —- 20 
4s Ca 76 26 
4s Se 64 — 
APPENDIX 


Numerical methods used to solve the equations (14) and (15) on BESK 


The Runge-Kutta method with constant intervals 6 was used. The equations 
were integrated inward, starting at a large radius R with a value P’(R). (Here 
the index of P;_,,.; is omitted for convenience.) Integrations were carried out 
with different starting values P’(R). The correct starting value P’ (R) is that for 
which the resulting P’(r) is zero at r=0. 

The following method was developed for the choice of P’(R): 

The first integration is carried out with a rough estimate of P’(R)=P(R), 
which value is to be given as a BESK parameter. The resulting function is de- 
noted Pj(r) with the value Pj (0) at r=0. The following flow-diagram will show 
how the different values of P’(R) are then automatically chosen. 

The procedure is continued until we get a function P),(r) for which P; (0)/ 
P), max is less than a certain number and until P;,,,(r) and P;, (r) are sufficiently 
close and have different signs for r=0. In Table 6, the method is demonstrated 
for the actual case, (2s—>p) for Li. 


Table 6. Choice of P’(R) for the actual case, (2s—>p) for Li. 


no. P’(R) sign of P’(0) 
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From be 
integer. Box A Orders are modified Integr. 
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the sign of Pj(O) is 


is performed in Box B 


tested 


To 
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ni+te1 niet 


cea. s 


p’ (R)=P’ 
mete 0% 


Fig. 1. Flow-diagram showing how the different values for P’(R) are automatically chosen. 
A is to be given as a BESK parameter. A is positive for the equations ( 14) and negative for 
the equations (15). 


The integral J = [Pr Pdr was calculated by Simpson’s formula for P;, and Pe 
0 
For the equations that were actually calculated, the value aunts was less 
n 


than 0.5%. 


To estimate how the choice of the constant interval 6 may have influenced 
the resulting value of I, I was calculated both for the constant interval 6 and 
Is —Is52 


for the constant interval It was found that I 
6 


was less than 0.5 % in 


the actual cases. 


Consequently the error in the values of « due to the numerical methods can 
be estimated to be less than 1%. 
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